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. . . Abstract 
(N ; 

^ ■ Multi-hadron production in inelastic neutrino-nucleon interactions is investigated 

^ . within the framework of the quark-gluon string model. The contributions of the 

O ! planar (one-Reggeon exchange) and cylindrical (one-Pomeron exchange) graphs to 

different observables is computed using a Monte Carlo program for the generation of 
hadrons produced from the decay of colorless quark- ant iquark strings. The suggested 
approach results in a satisfactory description of the experimental data on v{i')N — > 
li~{lJi^)hX reactions obtained recently at CERN by the NOMAD Collaboration. 
^ ! The data extends over a wide range of initial neutrino energies E^^ < 200 GeV/c 

I and momentum transfers 1 < Q < 7 GeV/c, well into the region where perturbative 

^ ' QCD calculations are not applicable. 

^ ■ PACS: 13.60Hb, 13.15.+g 

(N ■ 

■ The investigation of multi-hadron production in inelastic lepton-nucleon inter- 
^ • actions is a tool to study the dynamics of such processes and the nucleon structure 

ji^! function, particularly the dependence over a wide interval of values. In ana- 

Qj \ lyzing these types of reactions, the fragmentation properties of the struck quark in 

' deep inelastic scattering (DIS) in u^N and z/^A^ charged current events are often 

■ compared to those of the quarks produced in e~^e~ annihilation and ep inelastic 
^ ■ scattering. This comparison tests the universality of the quark fragmentation pro- 
cess in regions where non-perturbative QCD effects are important. Such an analysis 
has been performed recently in Ref. |jl[. It has also been found that a perturbative 
QCD calculation results in a good description of the charged hadron multiplicity 
< rich > measured in e~^e~ annihilation and ep DIS at initial energies E = ^/s = Q 
> 5 GeV /c [|T| , whereas a similar calculation of < rich > in the current region of pro- 
cesses z/(z/) p — » h X cannot reproduce the NOMAD data at Q < 5 GeV/c 
. This discrepancy is an important motivation for the present analysis of inelastic 
neutrino-nucleon interactions within the framework of the quark-gluon string model 
(QGSM) PI, ^, which is essentially equivalent to the dual parton model (DPM) de- 
veloped in Refs. 0, ^]. We calculate the different observables in a manner similar to 
Ref. 0, which focused on the pion multiplicity in z/(z/) A /i~(yU+) h X reactions. 
The main difference from the previous study is that we now use the Monte Carlo 
(MC) generation of hadrons produced from the decay of colorless strings. This pro- 
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cedure allows us to analyze all the observables measured experimentally and make 
a comparison with recent data 

In the QGSM p|, ^, hadron production in the reactions z/(z/) p h X \s 

described in terms of planar and cylindrical graphs, as shown in Fig. 1. The planar 
graph of Fig. 1(a) describes neutrino scattering off a valence quark, corresponding 
to the one Reggeon exchange in the t-channel P], whereas the cylindrical graph of 
Fig. 1(b) describes neutrino scattering off sea quarks in the proton corresponding to 
the one-Pomeron exchange in the t-channel @. The figures also show the occurrence 
of hadronization in the colorless quark-antiquark and quark-diquark strings. 

The relativistic invariant distribution of hadrons produced in the process z/(z/) p 
h X is defined as 

dN 

Pui.) H X - Eh^^^^^j^ , (1) 

where Eh and are the total energy and three momentum of the produced hadron, 
respectively; E' and Q are the energy and the solid angle of the final muon. It can 
be written in the following general form 

Pu(u) p^t,-(i,+) h X = ^{Q^) {Fp{x,Q^;z,pht) + Fc{x,Q^;z,pht)} (2) 

with 

$(Q2) = mE^' "^'w 



vr Q'^ + m 
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where G is the Fermi weak coupling constant, E is the energy of the incoming 
neutrino, m and mw are the nucleon and the ly-boson masses respectively, x = 
Q'^/2{pu ■ k) is the Bjorken variable, p^ and k are the four momenta of the initial 
neutrino and nucleon. Here Fp{x, z,pht) and Fc{x, Q^; z,pht) are the probability 
distributions of produced hadrons corresponding to the planar and cylindrical graphs 
of Fig. 1 respectively. The contributions Fp and Fc are computed by using the MC 
generation of all the quark-antiquark and quark-diquark strings drawn in Fig. 1(a) 
and 1(b). Generally, z is the light cone variable defined as z = {Eh +phz)/{E + Pz), 
where phz is the hadron longitudinal momentum respect to the initial neutrino having 
momentum p^, see fro example Refs. p, ^ and @, 3 . At large energies of final hadrons 
it is the longitudinal momentum fraction of hadron with respect to the neutrino in 
the rest frame of the proton-target. The variable z can be treated also as the 
Feynman variable xp = 2p*i^/Wx defined as the longitudinal momentum fraction in 
the hadronic center of mass system (h.c.m.s.). Here p2 is the longitudinal hadron 
momentum in the h.c.m.s., Wx is the mass of hadrons produced in the reaction, 
then pht is the transverse momentum of produced hadron with respect to the current 
(hadronic jet) direction, see for example Ref. [0]. 

The main ingredients for the calculations of observables in the discussed reac- 
tion are the quark distributions in a nucleon and their fragmentation functions to 
hadrons. In addition to the dependence upon x, quark distributions also depend on 
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and the transverse momentum kt- Following Ref. 0, we take a factorized form 
for these distributions: 



qf{x,Q'^;kt) = qf{x,Q^)gq{kt) , (3) 
with the function Qg chosen in the form 

9,ikt) = ^e-^'^ , (4) 

where B = 1/ < kt >~ A{GeV/c)~^, and < kt >~ 0.25GeV/c is the average 
transverse momentum of a quark in a nucleon. As for the function qf{x, Q^) we use 
the fit suggested in Ref. [^] including the true Regge x-asymptotic at x — >• 0, x — >■ 1 
and small Q^, and its QCD prediction at large Q^. More specifically, the valence u 
and d quark distributions according to Ref. [Q, have the following forms: 

/ n2 \ "^^^^ 



q,{x, Q') = - x)"(«^)+^ (^^^ j , (6) 

where 



We take the values of the constants from Ref. [Q: aji{0) = 0.4150 is the Reggeon 
intercept, 5„ = 1.2064, Ba = 0.1798, b = 0.6452, c = 3.5489 . 
The sea quark distribution in the proton is taken to be |§]: 

g_(x, Q') = Ax-^(«^)-^(l - (^^^ j , (7) 

where 

and A = 0.1502, a = 0.2631, Aq = 0.07684, Ai = 2.0, As = 1.1170. Note, this fit of 
the quark distributions in the proton describes all the experimental data from very 
small X up to X ~ 0.9 p. 

Generally, the fragmentation functions (FF) of quarks (diquarks) -D^j-^g^ into 
hadrons depend on the hadron momentum fraction Zi and the hadron transverse 
momentum p^t with respect to a quark (diquark) momentum direction. Here also 
we choose the factorized form 

Dg(gg){z,,Pht) = ( ^1 ) (^/.t ) , (8) 
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with the function gg again chosen in the form (Q). The functions D'^i^^^){zi) are 
constructed, according to the recursive cascade model procedure suggested in 0. 
They are found from the following integral equation 0: 



The function f{x) is chosen in the form 

fix) = x%i-xy (10) 

According to the main points of the QGSM, the fragmentation functions D^f^^^^^zi) 
should satisfy true Regge asymptotics at — 1 and zi — >• This limitation 
allows us to find the values of the parameters. The detailed procedure is presented 



in Ref. |]TOl- In general case the FF have to depend not only on zi,pht but also on 
Q^. At small they have to reproduce the true Regge asymptotic and at large 

they have to describe the e^e~ annihilation data. Since we analyze inelastic 
neutrino- nucleon interactions mainly at not large Q^, one can assume that the FF 
depends on it weakly and neglect this behaviour of D'^i^qq)- 

The MC computation results of the different observables in the reaction v + N ^ 
fi~ + h + X are presented in Figs. 2 to 5. The mean charged multiplicity in the 
current region is presented in Fig. 2. As it shown in [Q the NOMAD multiplicity 



values rich are very close to < rich > /2 results from e^e experiment |jTT| at energy 



E = y/s and < rich > from ep and i/^p at E = Q. In Fig. 2 we compare our results 
with the QCD calculation of the charged multiplicity from e^e~ < ri^^^ > /2. The 



open circles in Fig. 2 correspond to the QCD calculations ||12| , jTSl for the evolution 
of partons in the leading log approximation which give for < ri^^ > the following 
fit, see also [T^ , 



ngf^ = a + b exp{c^ln{QyQl)), (11) 

where a = 2.257, b = 0.094, c = 1.775, Qo = IGeV/c. The solid, long dash and short 
dash lines in Fig. 2 correspond to our Monte Carlo calculation at the initial neutrino 
energies E^ = 150GeV, E^ = A5GeV, E^ = 23.QGeV respectively. The NOMAD 
experimental data are averaged over the initial energy |jll|. It is seen from Fig. 1 
that the QCD fit can't reproduce the NOMAD data on the charged multiplicity at 
Q < 5GeV/c, whereas at large Q it describes these ones and < rich > from e~^e~, 
ep and u^p experiments satisfactory as it is shown in |1[]. The suggested approach 
describes the NOMAD data at Q < 5GeV/c rather well. 

In Fig. 3 the dependence of the inclusive muon spectrum in the process 
z/ + p — » + X is presented. The long and the short dash lines correspond to the 
the contributions of the cylindrical (one-Pomeron exchange) and the planar (one- 
Reggeon exchange) graphs to this spectrum, the solid line shows their incoherent 
sum. It is seen from Fig. 3 that the contribution of the one-Pomeron exchange 
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graph is decreasing with increase of Q"^ much faster than the one of the one-Reggeon 
diagram which dominates at > 20{GeV/c)'^ . 

In Fig. 4 the xp distributions of strange i^°-mesons and A-baryons are presented. 
The sohd hues in Fig. 3 correspond to our calculation, the experimental points 
are the NOMAD data 0]. This figure shows rather satisfactory description of the 
NOMAD data. 

In Fig. 5 our calculation results of the mean multiplicity of the strange hadrons 
Kg, A and A as the function of the neutrino energy with the NOMAD data 
are presented. 

In this Letter we have analyzed the multiple hadron production in inelastic 
neutrino-nucleon processes, The main conclusions of our work can be summarized 
as follows. By the analysis of the charged hadron multiplicity in the current region 
of the process h' + p^fi~ + h + Xa.tQ'^< 10{GeV/c)^ the nonperturbative QCD 
effects become very important. The conventional perturbative QCD calculation of 
< rich > doesn't reproduce the NOMAD data in this region, whereas the application 
of the QGSM or the DPM allows to describe these data rather satisfactory. The 
main contribution to this observable and the muon distribution 1 / Nda / dQ^ at small 

is coming from the cylindrical (one-Pomeron exchange) graph. Fig. 1(b), whereas 
at large the planar (one-Reggeon) graph. Fig. 1(a), dominates. Note, these two 
observables presented in Figs. 2,3 are related to the current fragmentation region. 
Recently in it has been shown that in the target fragmentation region these two 
contributions to the multiplicity of pions produced in the backward hemi-sphere in 
inelastic neutrino-nucleus collisions have the dependences opposite to the ones 
obtained in this paper. There is some analogous between the last reaction and the 
soft h — N processes at x^? ^ — 1 when the one-Pomeron exchange graph dominates 
at large initial energy, see for example 0. Therefore the contributions of the Fig. 
1(a) and Fig. 1(b) graphs have different dependences in the different kinematical 
regions, the current and the target fragmentation ones. 

The NOMAD data on the xp distributions of strange hadrons are described by 
the QGSM rather satisfactorily. It can mean that the fragmentation functions of 
quarks (diquarks) used by description of soft h — N processes are rather good to 
analyse the multiple hadron production in inelastic v — N interactions at not large 
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Figure 1: (a) The planar (one-Reggeon exchange) graph and (b) the cyhndrical 
(one-Pomeron exchange) graph. 




Figure 2: Mean multiphcity of charged hadrons in the current fragmentation region 
as a function of the momentum transfer Q. The open circles correspond to the 
QCD fit given by eq. (|TT|) ; the solid, long dash and short dash lines correspond to 
our calculations at = ISO.Gef , 45.Gef . and E^, = 23.6GeV respectively. The 
experimental points are the NOMAD data |l| 
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Q' (GeV/c)' 

Figure 3: The Q^-distribution 1/NdN/dQ^ of muons produced in the inclusive pro- 
cess u + p — > fi~ + X. The long dash and short dash lines correspond to the 
contributions of the cylindrical and planar graphs respectively. The solid line is the 
sum of these contributions. 
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Figure 4: The xp distribution of strange hadrosn. The upper panel corresponds to 
the distribution dN/xp of mesons in semi-inclusive process v+N — *■ ix~ +K^+X, 
whereas the down panel shows the xp distribution of A'' hyperons produced in 
V + N /x^ + A° + X reaction. The experimental points are the NOMAD data . 
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Figure 5: The average multiplicity of the strange hadrons as a function of the 
neutrino energy Ei^. The upper line corresponds to < n > of the mesons, the 
middle line is the mean multiplicity of A° hyperons and the down curve shows < n > 
of the antilambda hyperons A°. 
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